Quasiparticle spectrum inside the vortex core: 
crossover from dirty to clean limit 
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The quasiparticle spectrum inside the vortex core in the mixed state of a strongly type-II super- 
conductors are studied. The s-wave symmetry for the gap parameter is assumed. The crossover 
behavior from dirty to clean limit is shown by numerical calculation based on the random matrix 
theory. 
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The study of electric structure of vortices is very im- 
portant for understanding the behavior of the supercon- 
ductors in a magnetic field. The quasiparticle energy 
levels inside the vortex core is firstly studied by Caroli, 
de-Gennes, Matricon jj]]. In the two dimensional case, 
the discretized energy levels with equal spacings are ob- 
tained and the quasiparticles are localized in these lev- 
els. The existence of these localized states is confirmed 
experimentally on NbSe2 f|||- The scattering between 
the levels due to impurities inside the vortex broadens 
the discretized energy levels and the system becomes ex- 
tended. In the clean case, where the number of impurities 
are moderately small, the periodic spectrum is obtained 
. In the dirty limit the level statistics inside the vor- 
tex core is studied using the random matrix theory [MB . 
In the latter limit, the impurity potential is given by 
the random symplectic matrix and the function form for 
the density of states near the Fermi energy is obtained. 
The crossover behavior between these two limits is re- 
cently investigated analytically with the random matrix 
theory [pj where the impurity potential is handled by a 
random symplectic matrix coupled to an external non- 
random source matrix |hJ. 

In the extremely type-II layered superconductors, the 
vortex core is treated as a normal disk of the radius £o 
inside a superconductor where £o is the coherence length. 
In the polar coordinate, taking the adequate gauge con- 
dition, the angular dependencies can be eliminated and 
the Bogolubov-deGennes equation is solved with the ba- 
sis which is written by the Bessel function (r) . Due to 
the s-wave symmetry of the conventional superconduc- 
tors, the eigen modes are fully determined only by the 
index of the angular momentum denoted by This one- 
dimensional feature of the quasiparticle modes in a vortex 
makes the numerical study for the level statistics inside 
the core very simple. When the short-ranged impurity 
scattering is introduced, the Hamiltonian is given by an 
infinite matrix. The size of this matrix corresponds to 
the number of eigen modes labeled by /i (— oo < \l < oo) 
and the matrix elements are given by the explicit wave 
functions for the unperturbed eigenstates. The density of 
states is obtained by numerical calculation of the eigen- 
values for an NxN scattering Hamiltonian at large value 



of N introducing some adequate cut-off for the eigen en- 
ergies. 

In this work, we present the numerical calculation for 
the crossover behavior of the density of states (DOS) in- 
side a vortex core from the dirty to the clean limit in 
the presence of impurity scattering between the levels. 
The two limiting cases, superclean and dirty case, are 
determined by the scattering time r between the Caroli- 
deGennes-Matricon (CDM) excitation states inside the 
vortex core. In the superclean case coq 3> 1/t where ujq 
is the level spacing of order u>q ~ A 2 /Ep. In the dirty 
case 1/t » A. As shown by Koulakov and Larkin [Q, 
the spectrum becomes an luq periodic function of energy 
in the wide region of the intermediately clean case. We 
define the dirty limit where the number of impurities are 
fairly large 3> 1 and the random matrix approach for 
the impurity scattering is appropriately applied p0[ . 

The Bogolubov-deGennes equations for the two com- 
ponent excitation wave functions ip — ( u , v ) inside the 
vortex of the s-wave superconductor is written as 



2m dr 2 



E F + V(r))+a x ReA(r) 



<T y lmA(r) 



ijj = Ei\) 



(1) 



where a x ,a y ,a z are Pauli matrices and A(r) is the or- 
der parameter. V(r) is the disorder potential which is 
produced by the short-range impurities 



V(r)=J> i <i(r-r i ) 



(2) 



where Vi, rj are the strength and the position of i-th 
impurity, respectively. The summation is taken over all 
impurities. We assume that the magnetic field is weak 
(B <C H C 2) and the order parameter is given by 



A(r) = A(r)e* 



(3) 



where r = |r| and = Arg(r). The CDM excitation 
spectrum is obtained for V(r) = case as 



(4) 
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where C is the normalization constant, K(r) = Ar/up = 
r/^n and // = 0, ±1, ±2, • • ■. The excitation energies are 
equidistant as 



(5) 



(We take h = 1 hereafter.) We assume that the gap A is 
large enough and there are a large number of excitation 
levels N ~ A/cj 3> 1- When the scattering between the 
CDM levels due to the impurity potential is present, the 
wave function is given by a linear combination of each 
levels 



(6) 



and the excitation spectrum E is obtained by following 
equations 

det[diag(£0-£) + ]Ti ! ] =0 (7) 

A* nm = Ce^(™-")[J„(fc F r l )J m (fcf n) 

-J n -i{kFri)J m -x{k F ri)] (8) 

where = (rj,#i) is the position of the i-th impurity. 
Since k F ri 3> 1 in our clean limit, we can use an asymp- 
totic expansion for the Bessel functions in the matrix 
elements A* , 



Q e i(m-n)e, gm [ 2k F r. 



(9) 



The density of states is obtained by solving Eq. (0) nu- 
merically and averaging over the positions of impurities 
which are randomly generated. Fig.l shows the result 
for the case where the number of impurities Nj = 4, the 
number of CDM basis N = 50 and the averaging has been 
done over 6000 generations of the impurity positions. 
The particle-hole symmetry of the spectrum is evident 
as seen from the plot, and it shows ujq periodic structure 
which can be seen from the enlarged plot near the Fermi 
level (E — 0). The analytical result by Koulakov and 
Larkin is written as 



2 9 (kE 
P(E) = — sin 2 — 



(10) 



Our result of Fig. 1(d) coincides with this analytical re- 
sult. We obtain relatively sharp peaks when the matrix 
elements are given by the Bessel functions as Eq.(||). 

In the dirty limit, it is obvious that the scattering 
Hamiltonian becomes an complex symplectic random 
matrix from the explicit form of A l nm in Eq.(^). We take 
2iV-CDM basis near the Fermi level (-N + 1 < fi < N) 
and rearrange the order of raws and columns of the ma- 
trix A in order to divide A into four blocks , in each of 
which the indices are even-even, even-odd, odd-even and 
odd-odd, respectively. Since J- n {z) = (— l) n J n (z), we 
have 
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FIG. 1. The density of states near the Fermi level for the 
clean limit. The number of CDM modes is 50 and the number 
of impurities Ni = 4. The eigenvalues of the scattering ma- 
trices is evaluated with the FORTRAN package "EISPACK" 
and impurity averaging has been done over 6000 iterations. 
For the figures (a),(b) the matrix elements are given by Eq.(8), 
and for (c),(d) it is given by Eq.(9), respectively. 



A = 



(11) 



where a is an N x N hermitian matrix and b is an N x N 
symmetric complex matrix. The symplectic structure of 
A is evident from the relation 

A T J + JA = (12) 

where J is written with the N x TV identity matrix / 



J 



I 

-I 



(13) 



In Fig. 2 the DOS for the dirty limit is shown which is 
evaluated by the calculation of the eigenvalues of the 
Gaussian symplectic random matrices of size 100. It 
is consistent with the conjecture of analytical study of 
the level statistics for class C of the random matrix the- 
ory. When the Hamiltonian is given by a random matrix 
taken from a symplectic ensemble Sp(N), the DOS near 
the Fermi level is given by 



P(E) = 1 - 



sin(27rS/a;o) 
2ttE/uj 



(14) 
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FIG. 2. The density of states near the Fermi level for the 
dirty limit. The number of CDM modes is 50. (a) The full 
range plot, (b) The enlarged plot near the Fermi level. The 
solid line is the analytical result Eq. (14). 

This result is also shown by the solid line in Fig. 2(b). 
Although the random matrix theory is fully phenomel- 
ogical and its applicability to the real system is left un- 
determined, the investigation of the crossover behavior 
between the above two limits by means of the matrix 
theory is quite important since it may provide a lot of 
interesting behavior of the level statistics inside the vor- 
tex core which can be measured by the STM experiments 
directly. This level statistics also determines the trans- 
port properties of superconductors in a magnetic field. 

In order to investigate the crossover behavior of the 
spectrum between the above two limiting cases, we in- 
troduce the effective Hamiltonian in the form 



TL = Tin 



V 



(15) 



where TLq is the Hamiltonian for the clean limit whose 
eigen energies arc written as 



(16) 



where z is some function of the position of impurities. 
The spectrum is obtained by averaging over z. V is 
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FIG. 3. The density of states near the Fermi level for var- 
ious value of R. The number of CDM modes is N = 50. 
The energy is scaled by both 7? and u>o- The averaging over 
impurities has been done for 6000 iterations. 



the random symplectic matrix whose probability distri- 
bution is assumed to be Gaussian. We also introduce a 
parameter R which determines the strength of the exter- 
nal source Tio respective to the random potential V . The 
matrix elements of TCq is written by 



Ni 



Re ie ^ m -^[J n (k F n)J m (k F n) 

-J n -i{kFri)J m -i{k F ri)] 



(17) 



In the numerical calculations, the rearrangement of the 
indices of the matrix element A l nm which generate the 
four blocks, even-even, even-odd, odd-even and odd-odd 
is applied. We fix the amplitude of the random matrix 
in the calculation. When R — > oo, the clean limit is 
obtained and taking R = gives the dirty limit. We 
obtain the DOS for several values of R in Fig. 3. The 
vertical axis is the scaled energy E/(Ru>o). We see that 
the peaks grow as the value of R is increased. 

In Ref. p0| , the crossover from the clean limit to the 
dirty limit is investigated with the case where the deter- 
ministic matrix with the eigenvalues Eq.([l^) is coupled 
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to the random matrix. By means of the generalization of 
the unitary ensemble matrix model, they discussed the 
symplectic ensemble case. The explicit formulae for the 
density of states is given by 



P{E) 



N 
C 7 
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_ 4N 2 /C 2 ANE, 



cos(-^-) 



A(E 2 + N 2 /C 2 ) 



e~ 4N ' c sin( — ) (18) 



where the parameter C determines the strength of the 
external source Tio respective to the random potential V 
and N is the matrix size. Our numerical result is qualita- 
tively in good agreement with this analytical result near 
the Fermi level. 

On the other hand, the crossover behavior can be ob- 
served by increasing the number of impurities Ni directly. 
In this case, we take the Hamiltonian as 



E° 



N 



O 



O 



E° 

N 



(19) 



where A 1 is the scattering matrix Eq.(g) generated by 
one impurity at site = (r*j,0j). Fig. 4 shows the re- 
sult of the numerical calculation for various values of N . 
As seen from Fig. 4, the density of states near the Fermi 
level becomes non-vanishing except at the origin when 
the number of impurities is increased. The spectrum near 
the Fermi energy for the large value of Ni is similar to 
the dirty limit case where the Hamiltonian is given by 
the symplectic Gaussian random matrix. This crossover 
behavior appears only in the narrow region of the energy 
levels around the Fermi level. 

To summarize, we have studied the quasiparticle spec- 
trum inside the vortex core numerically for the clean, 
dirty and crossover case of the conventional (s-wave) su- 
perconductors in a magnetic field. By means of the ran- 
dom matrix theory, the crossover behavior of the vor- 
tex density of states is obtained. In the relatively dirty 
crossover region, the density of states becomes zero only 
at the Fermi level and the spectra near the Fermi level is 
obtained numerically which is consistent with the analyt- 
ical result. The spectrum for the dirty limit is reproduced 
by increasing the number of impurities in the scattering 
matrix in the energy region near the Fermi level. The 
investigation of the vortex quasiparticle spectrum in the 
case of the unconventional superconductors such as d- 
wave superconductors would be very interesting which is 
in preparation. 

The author acknowledges helpful discussions with Shi- 
nobu Hikami. The numerical calculation were made by 
using the Vector Parallel Processor, Fujitsu VPP300 at 
CCSE. 
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FIG. 4. The crossover behavior of the density of states near 
the Fermi level obtained by increasing the number of impu- 
rities Ni from the super clean limit Ni — 1 to the dirty case 
N = 100. The number of CDM modes is N = 50. 
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